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High Frequency electron paramagnetic resonance has been used to observe the magnetic dipole, A 
Ms = ± 1, transitions in the S = 9 excited state of the single molecule magnet FegBrg. A Boltzmann 
analysis of the measured intensities locates it at 24 ± 2 K above the S = 10 ground state, while the 
line positions yield its mag netic parameters D = -0.27 K, E = ±0.05 K, and B2 = -1.3 x 10"*^ K. 
D is thus smaller by 8% and E larger by 7% than for S — 10. The anisotropy barrier for S = 9 is 
estimated as 22 K,which is 25% smaller than that for 5 = 10 (29 K). These data also help assign the 
spin exchange constants(J's) and thus provide a basis for improved electronic structure calculations 
of FesBrg. 

PACS numbers: 75.50.Xx, 75.60.Jk, 75.75.+a, 76.30.-v 



I. INTRODUCTION 



Single molecule magnets (SMM's), defined as com- 
pounds where a magnetic domain can, in principle, 
be reduced to a single moleculei*^, have recently been 
of high theoretical and experimental interest due to 
their novel properties and potential applications. These 
properties and potential applications include quan- 
tum tunneling of their magnetization (QTM)^i^i^, 
whose detailed mechanism is still not fully understood, 
molecular memory devices^iS and elements of quantum 
computers^. One of the best characterized SMM's 
is [(C6Hi5N3)6Fe8(/i3-0)2(/^2-OH)i2]Br7(H20)Br-8H20, 
abbreviated FegBr^Siii, whose main spin-bearing skele- 
ton is shown in Fig. ^ 

Studies by magnetization^^, neutron 
scatteringi2ii4, electron paramagnetic resonance 
(EPR)"d5,i6ii7,i8ii9Wi^ NMRS^ techniques 

have established that the ground state has a spin value 
S = 10. Figure n shows a schematic of the spin config- 
uration of the Fe'^+ (5 = 5/2) ions. At temperatures 
below 1 K, the magnetization relaxation takes places 
via QTM^. While a great deal of progress has been 
made in understanding the nature of QTM in FesBrg, 
many questions still remain unclear. For example, the 
magnitudes of calculated tunneling rates are much lower 
than the observed values^. Secondly, there is a lack of 
data on the nature and magnitude of the spin-exchange 
constants, J's, between the eight Fe'^+ ions in the FegBrg 
core. The best estimates come from the temperature 
dependence of the dc magnetic susceptibility, Xd<~^- 
The Xd6^^ could be fitted by several sets of J's. One 
of the criteria of such a procedure is the prediction of 
the proper spin S value of the ground state, together 
with the location of the excited states, in particular 
the 5 = 9 state. The Xdc ht yielded at least two sets 
of J's, but the set providing the better fit yielded the 
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FIG. 1: Schematic representation of FegBrg— . The arrows 
represent spin orientations of the Fe^+ ions in the 5* = 10 
ground state. The organic ligands and the eight Br~ anions 
have been omitted for clarity. Each Fe'^"'' has S = 5/2, thus 
the ground state can be seen asS = 6x5/2 — 2x5/2 = 10. 



position of the S = 9 state to be less than 0.5 cm~^ 
above the ground stateiS. The other set predicted the 
5 = 9 state at greater than 25 cm^^ (36 K) above 
the ground state. Clearly, from the point of view of 
understanding the magnetic structure of FcgBrg, it is 
important to experimentally determine the location of 
the excited states, in particular the 5 = 9 manifold. 
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A peak around 30 cm"'^ has been reported in the far 
infrared spectrum^^, but it was not estabhshed whether 
it was a normal molecular vibration involving metal 
ions, or a transition to the S — 9 level from the 5* = 10 
ground state. Moreover this peak was not observed for 
the closely similar compound Fe8Br4(C104)4, whose spin 
Hamiltonian (Eq. parameters are known to be close 
to those of FcsBrs^l 

The magnetic properties of the ground state of FcgBrg 
have been well described by the following spin Hamilto- 

nian, with S = 10ll,15,16,17,18,19,20,21,25,26,27. 

H = fiBB-g-S+DSl+E{S^^^Sy)+Bl04+Bi0i+Bt0t, (1) 



where the first term is the Zeeman interaction, D rep- 
resents the usual zero-field, uniaxial anisotropy param- 
eter, and E the second order rhombic anisotropy. The 
fourth order terms are given by: O4 — 35 - [30 
SiS+iy 25]Sl - 6S{S +l) + 3 S'^iS + if, Oj 1/4[7S'2 
- S{S+1) - 5]{Sl + 52 )+ (^2 ^ 52 )[752 _ ^(^^1) _ 

and Of = 1/2(54 + St). 

Here, we report on an EPR detection of the S = 9 state 
of FcgBrg, and preliminary evaluation of its spin Hamil- 
tonian parameters as defined in Eq. We observe 
a series of peaks, called /3 transitions, which we assign 
to an 5 = 9 excited state of FegBrg. From the temper- 
ature variation of the line intensities, we establish that 
the 5 = 9 manifold lies at 24 ± 2 K above the S* = 10 
ground state, in contrast to the suggestion of > 36 K 
from susceptibility analysisi^. 

Section II below describes the single crystal prepara- 
tion and the EPR instrumentation. The results obtained, 
and their analysis, are presented in section III, with the 
discussion presented in section IV, and the main conclu- 
sions presented in section V. 

II. EXPERIMENTAL 

[(C6Hi5N3)6Fe8(/Z3-0)2(/i2-OH)i2]Br7(H20)Br-8H20 
(FegBrg) was synthesized following the procedure of 
Weighard^. Relatively large (2 mm x 2 mm x 0.5 mm), 
optical quality single crystals were prepared by slow 
evaporation. The crystals were aligned with the Zeeman 
field applied along the easy axis of magnetization by 
sight, to within a few degrees^^. The final orientation 
was confirmed by EPR splittings, being the extremum 
for the canonical orientations. The EPR measurements 
were made using a variable frequency (44 - 200 GHz), 
cavity based, high sensitivity spectrometer described 
earlierii^iiS*2&. The main component of the spectrometer 
is a millimeterwave vector network analyzer (MVNA), 
a phase-sensitive, fully sweepable, superheterodyne 
source-detection system. A variable fiow cryostat situ- 
ated within the bore of an 17 T superconducting solenoid 
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FIG. 2; EPR spectrum of FegBrg at 131 GHz and 35 K 
with the zeeman field applied along the easy axis (bottom 
panel), along with the energy level diagram corresponding to 
the S = 10 spin system (top panel). The S = 10 energy level 
diagram has been constructed through the spin Hamiltonian 
parameters of CaciufFo et alw—. The single-headed arrow at 
0.8 tesla designates the expected position where a /3io would 
appear if the /3 transitions originated from an 5* = 10 state. 



allows for temperatures down to 1.5 K, with an accuracy 
of ± 0.01 K. The high sensitivity of the MVNA technique 
(10^ spins G^^ s~^) allows for observation of the low 
level transitions of the FegBrg ground and excited states 
in a single crystal, and for angular variation studies, as 
described earlieriSiiSiSi. 



III. RESULTS 

Figure [3 (bottom panel) shows a typical EPR spec- 
trum of FegBrg at 131 GHz with the Zeeman field applied 
along the easy axis of an FegBrg single crystal at 35 K. 
The spectrum consists of a series of strong peaks a_io, 
a_g, a_g, etc; the subscripts represent the spin projec- 
tion quantum number, Mc,, corresponding to the level 
from which the EPR absorption transition originates, in 
the S* = 10 ground state following the convention intro- 
duced earlieriiiSiiZ.. The a transitions have been very 
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well analyzed earlieniiiiSiii, and have been shown to arise 
from the magnetic dipole (AM^ = ± 1) transitions within 
the 21 Ms levels of the S* = 10 multiplet. The specific 
peak assignment is indicated in the top panel for the a 
transitions. In addition to the a resonances, there are ad- 
ditional peaks present in Fig. |21that are labelled as the /3 
transitions and are the focus of the current investigation. 

Our analysis procedure consisted of three steps: (a) 
to ascertain that the /3 transitions are from an excited 
state, (b) to determine the spin multiplicity of this ex- 
cited state, and (c) to deduce the spin Hamiltonian pa- 
rameters for the P spin system and its energy position 
relative to the ground state {S — 10). 

Direct evidence that the (3 transitions originate from 
a thermally populated excited state is provided by the 
temperature dependence of their intensities. Figure |31 
shows spectra at 5, 15, and 35 K, respectively, giving 
clear experimental evidence that the /3 transitions arise 
from a thermally populated excited state, because their 
intensities rapidly decrease as the temperature is lowered. 

In order to quantitatively measure the intensities of 
the /3 peaks as a function of temperature, the spectra of 
the ground and excited states needed to be separated. 
The separation was accomplished by using a Gaussian 
fit for each individual peak. The validity of Gaussian 
fits, especially for the low field transitions, has been pre- 
viously established^^. The criterion for the goodness of 
the fit was that the sum of the a and /? transitions mir- 
rors the experimental data quite well by visual inspection 
and also by minimizing the remaining intensity obtained 
by subtracting the Gaussian fits from the experimental 
spectra. This remaining intensity was within the experi- 
mental noise in our separation procedure. Each Gaussian 
fit was then summed to yield a separated 5 = 10 spec- 
trum. The spectra were then separated by subtracting 
the Gaussian fits of the a peaks from the experimental 
spectra. These Gaussian fits were then summed, giving a 
separated S* = 10 spectrum. The remaining (3 peaks were 
then fit with Gaussian functions, and summed, yielding 
a separated S = 9 spectrum. We verified that these sep- 
arated spectra agreed well with the experimental data by 
taking the sum of the a and /? spectra and comparing it 
with experiment, as shown in Fig. 0Ja). 

Figure 01 also shows the relative decrease in the inten- 
sity of the separated excited state (/?) spectra in relation 
to the 5* = 10 spectra as the temperature is decreased. 
The spectral envelope present in the S* = 10 spectra at 
5, 15, and 35 K is generally evident in the corresponding 
excited state spectra, though the relative intensities are 
modified slightly due to differing matrix elements in the 
transition probabilities, Pnm oc |< ■0ri I "5*+ | ipm >P, of 
the two spin systems, given by Eq. (jJJ: 

P(x[S x{S + l)- Ms X (Ms -I- 1)]. (2) 

In order to analyze the intensity of the (3 peaks we 
normalized their intensities. The normalization process 
involved dividing the intensity of a /? peak, with a given 
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FIG. 3: Temperature dependence of EPR spectra of FegBrg, 
for B II z at 131 GHz and 5, 15, and 35 K. 



Ms value, by the corresponding a peak with the same Ms 
value. This enabled us to ignore instrumental effects and 
spectral envelope changes with temperature. At 5 K, the 
(3 peaks are barely discernable from the noise level of the 
spectrum. At 10 K, however, they increase in intensity 
enough to emerge from the a transitions. The decreas- 
ing intensity of the (3 peaks, upon lowering temperature, 
unambiguously designate the (3 peaks as being due to a 
thermally populated excited state. 

Once it was concluded that the (3 peaks originated 
from an excited state, the spin multiplicity of the ex- 
cited state needed to be determined. The presence of 
an excited state (with perhaps S = %) close to the 
ground state of FesBrg has been previously inferred from 
magnetic susceptibilitjii^, and mentioned in subsequent 
^SR^*», EPRi2&, and neutron diffraction studieE22,, with- 
out any evidence for its location or multiplicity. Herein, 
the spin multiplicity of the excited state has been de- 
termined by two independent methods. First, by the 
location of the leading peak (/S-g) in the set of peaks 
assigned to the S" = 9 state. Figure |21 shows an experi- 
mental spectrum taken at 131 GHz and 35 K. As is evi- 
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FIG. 4: (a) Experimental, resummed a and /3, separated S = 
9, and separated S = 10 spectra at 35 K and 131 GHz for B 
I z. (b) Separated spectra at 15 K, (c) and at 5 K. 



dent in all spectra taken, there is no (3 peak between the 
a-io and a^g peaks. Furthermore, the a_io transition 
is symmetric and shows a clear Gaussian shape, as would 
be expected at this temperature for a well separated, in- 
dividual peak^^. This lack of a P^iq transition, barring 
spin Hamiltonian parameters being very different from 
those for the ground state, is strong evidence that the 
excited state is 5 = 9. The location of the first f3 peak 
is consistent with slightly modified spin Hamiltonian pa- 
rameters and a spin multiplicity of 5 = 9 for the excited 
state, as anticipated theoreticallj^. 

Additional, more quantitative support, that the spin 
S of this excited state is indeed S = 9 was provided by 
computer simulations which were run using SIM^*^. The 
procedure was first checked for the 5* = 10 state for which 
the parameters are knowni2iiL2^. The simulations for the 
S = 10 state were performed with the spin Hamiltonian 
parameters previously determined by Caciuffo et ali^ us- 
ing neutron scattering. A typical comparison is shown in 
Fig. O The simulated spectra for the 5 = 10 ground 
state were in close agreement with our experimental re- 
sults, thereby validating the simulation procedure. 

The /3 transitions were then accurately simulated for 



TABLE I: Experimental and simulated peak positions at 131 
GHz for the 5 = 9 state, for B || z, using D = -0.27 K, E = 
±0.05 (±0.015) K, B4'' = -1.3 X lO"** K ,B4^ = K, B4'' = 
K 
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FIG. 5: Experimental and simulated EPR spectra of FesBrg, 
for B 11 z at 35 K and 131 GHz. 



the three frequencies utilized (110, 131, and 155 GHz), 
using the spin Hamiltonian parameters obtained in the 
present study. A field-frequency plot is shown in Fig. 
El Quite good agreement can be seen between the sim- 
ulated curves (continuous lines) and the observed peak 
positions. A more explicit and quantitative comparison 
is shown in Table I, using the measured data for 131 
GHz. The D parameter decreases in magnitude by 8%, 
to -0.27 K, while the E term increases by 7%, to ± - 
0.05 K (± 0.015 K). Surprisingly, the term is similar 
in magnitude, -1.3 x 10"^ K, but opposite in sign to 
the S" = 10 parameter (1.01 x 10"^ K). Due to the fact 
that the experimental spectra were taken along the easy 
axis of magnetization, the B4 and B4 terms were not 
included in the simulations. A simulated spectrum at 
35 K, with the spin Hamiltonian parameters determined 
for both the S' = 10 and 5" = 9 states is shown in Fig. 
|S1 Linewidths were not an optimized parameter in the 
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FIG. 6: Frequency dependence of the first eiglit /3 transitions 
in FesBrs with B || z. The solid lines portray the frequency 
dependence given by the spin Hamiltonian parameters deter- 
mined in this work, whereas the solid squares are the observed 
field positions. 

simulations, though the relative intensities matched the 
experimental data well. Therefore, the simulated spec- 
trum presented in this work is the sum of the S" = 10 and 
S — 9 separated spectra using peak positions generated 
by SIM^. The agreement between the experimental and 
simulated spectra, as shown in Fig. |31 can be seen to 
be quite satisfactory, thereby supporting the parameter 
assignment. We are thus able to assign the full experi- 
mental spectrum to transitions in the S = 10 multiplet 
(a's) and the S = 9 multiplet (J3's), as shown in Figure 

m 

Once the spin Hamiltonian parameters had been deter- 
mined, the relative intensities of the a and f3 transitions 
were used, at temperatures from 5 to 35 K, to determine 
the location of the S = 9 state above the ground state. 
The intensity of a specific peak between given Mg and 
Ms+i states, is proportional to the population difference 
between the Ms and Mg+i states and the transition prob- 
ability, P, as given in Eq. Q: 

/ cx P(iVM. - iVM.+ J/^. (3) 

Therefore, assuming very similar partition functions 
and D values for the two states, the intensity ratio be- 
tween two transitions of the same Ms states in the 5 = 9 
and 5 = 10 manifolds, respectively, is given by Eq. 

J9//10 = iP9/Pw)exp{~AEw^9/kT), (4) 

where k is the Boltzmann constant and AEio-g is the 
energy difference between given Ms states in the S = 
10 and S = 9 manifolds. The areas of the Gaussian 
fits, for a given Ms to Ms+i transition, were factored 
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FIG. 7: (a) Temperature Dependence of the intensity of the 
-8 to -7 transition in the 5 = 9 state, I9, normalized to the 
-8 to -7 transition in the 5* = 10 state, Iio. The curve joining 
the experimental points is a guide to the eye. (b) Boltzmann 
analysis of the normalized intensities of the 5* = 9 spin state. 
The slope yields the excitation energy as 18 ± 2 K between 
level with the same Ma values in the S = 10 and S = 9 
manifolds. Addition of the zero field splitting between the 
Ms = — 10 and Ms = —9 levels (5.5 K) leads to the location 
of the S — 9 state at 24 ± 2 K above the ground state. 



by their transition probabilities, in both the 5 = 9 and 
5 = 10 manifolds, in order to determine their ratios. 
Due to the fact that each /3 transition is normalized to its 
corresponding a transition, the ratio of intensities should 
be constant, regardless of the specific Ms pair, for any 
given temperature. 

A Boltzmann analysis of the intensity ratios is shown 
in Fig. [7fb). The ratios of intensities of the (3 to a 
transitions were compared at 5, 15, 20, 25, 30 and 35 K, 
and plotted vs. inverse temperature (T~^). The slope 
yields the energy difference (18 K) between a given Ms (9) 
state and the corresponding Ms (10) state. Therefore, 
the energy difference between the Ms = 10 and Ms = 
9 in zero field (5.46 K) must be added to the energy 
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obtained from the Boltzmann analysis, yielding an energy 
difference A of 24 ± 2 K (16.7 ± 1.5 cm-i). Figure El 
shows a schematic of the energy levels in zero field for the 
S" = 10 and 5 = 9 manifolds based on the present study. 
Clearly, the higher levels of the S = 9 state overlap 
with the lower Ms levels of the 5=10 state, indicating 
at least a partial breakdown of the single spin model. 



IV. DISCUSSION 

Our determination that the 5 = 9 manifold in FcsBrg 
is located 24 ± 2 K above the 5 = 10 ground state is 
in contrast to with earlier suggestions based on magnetic 
susceptibility (> 36 K)i^. The four central Fe'^+ ions 
in the FcsBrs core (Fei^+, Fe2^+, Fe3^+, and Fe4^+ (as 
shown in Fig. 1)) can be described as that of a butter- 
fly configuration. All magnetic coupling interactions, J, 
between Fe'^+ ions in FcgBrg have been determined to be 
antiferromagnetioi^. Therefore, the magnitude of these 
coupling constants dictates both the location and spin 
topology of the 5 = 9 excited state. 

As mentioned in the Introduction, a detailed analysis 
of the dc magnetic susceptibility data led Delfs et a/.— 
to two reasonable sets of exchange parameters; (a) J1-2 
= 20 cm^^, Ji_3 = 120 cm~^, Ji_5 = 15 cm^^, and J3_5 
= 35 cm-\ and (b) Ji^2 = 102 cm^^, Ji_3 = 120 cm"i, 
Ji_5 = 15 cm~-^, and J3_5 = 35 cm^^. While set (b) 
provided a much better fit to the experimental data, it 
predicted the position of the first excited state, 5 = 9, 
at less than 0.5 cm~^ above the 5 = 10 ground state. 
We do note, however, that Delfs et al. did not include 
any zero field splitting terms in their susceptibility anal- 
ysis. Nevertheless, this same basic configuration of ex- 
change constants has been recently supported by detailed 
symmetry-based calculations by Raghu et al.^K Though 
the magnitudes of the coupling constants calculated by 
these authors^^ arc different from those of DelfsA^., the 
dominance of the Ji_3 interaction over other magnetic 
couplings remains consistent. The present study sup- 
ports the essential correctness of set (a) . 

The coupling set (a) of Delfs et alm^ and the best set 
of Raghu et alM, both show that Ji_3 dominates the 
exchange interactions. The perturbation leading to the 
5 = 9 excited state must result from the smallest dif- 
ference in J's acting on the same ion or symmetrically 
equivalent set of Fe'^^ ions. Thus it seems reasonable 
to deduce that this perturbation does not involve Ji_3, 
hence the butterfly core (Fci, Fe2, Fe3, Fe4), but rather 
the Fe's on the corners of the cluster (Fe5, Fcg, Fey, Feg) 
via some linear combination of their wavefunctions. 

The spin Hamiltonian parameters determined for the 
5 = 9 state provide some insight to the origin of the 
anisotropy of the cluster. The 5 = 9 parameters are 
slightly different from those of the 5 = 10 ground state. 
The 7% larger E value for 5 = 9 is in accord with in- 
creased transverse distortion in the Fcg structure. The 
decrease in D with decreasing magnetic moment indi- 
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FIG. 8: Schematic for the energy levels of both the S = 10 
and the S = 9 states in zero magnetic field. S = 9 is located 
at an energy A = 24 ± 2 K as marked. 



cates that the anisotropy present in the FcgBrg core has 
some dipolar contribution rather than arising purely from 
a spin-orbit interaction. Similarly, the change in sign of 
the B4 term indicates that B4 originates from many-body 
interactions between Fe^+ ions, and not from a collective 
sum of individual B4 terms. This argument is in line with 
the fact that a B4 term needs an effective interaction in- 
volving at least four spins. Alternatively, this significant 
change in B4 may be a result of the breakdown of the 
single spin model, as has been proposed by Katsnelson et 
al. in connection with Mni2-Acetate2^. Additional, de- 
tailed angular variation studies are underway for precise 
measurement and understanding of these questions. 

The anisotropy barrier, estimated from D and E values 
for the 5 = 9 manifold is 22 K, as compared to that for 
the 5 = 10 ground state (29 K). 



V. CONCLUSIONS 

Using variable frequency, high-field, EPR measure- 
ments on single crystals of FegBrg, we have detected a 
set of transitions, labelled as f3i, which have been con- 
clusively assigned to the 5 = 9 spin multiplet, located 
at an energy A = 24 ± 2 K (17 ±1.5 cm~^ ) above 
the ground (5 = 10) state. The spin Hamiltonian pa- 
rameters have been determined to a good accuracy, and 
differ from those of the 5 = 10 state: D is smaller by 
8%, while E is larger by 7%. These parameters yield the 
anisotropy barrier (~ DS^^ ~ 22 K), about 25% smaller 
than for 5 = 10. 84° for 5 = 9 also shows a dramatic 
change, the sign is opposite to that for the 5 = 10 state. 
Although electronic structure calculations have been re- 
ported for FcgBrgiSiiSiM, there has been little definitive 
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data on excited states. The results of the present study support, 
should serve as sensitive basis for more refined theoreti- 
cal modelling of the bonding and magnetic properties of 
these materials. 
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